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ABSTRACT 



Context. Recent discoveries of intensity correlated downflows in the interior of a sunspot penumbra provide direct evidence for 
overturning convection, adding to earlier strong indications of convection from filament dynamics observed far from solar disk center, 
and supporting recent simulations of sunspots. 

Aims. Using spectropolarimetric observations obtained at a spatial resolution approaching C/l with the Swedish 1-m Solar Telescope 
(SST) and its spectropolarimeter CRISP, we investigate whether the convective downflows recently discovered in the Ci line at 
538.03 nm can also be detected in the wings of the Fe i line at 630.15 nm 

Methods. We make azimuthal fits of the measured LOS velocities in the core and wings of the 538 nm and 630 nm lines to disentangle 
the vertical and horizontal flows. To investigate how these depend on the continuum intensity, the azimuthal fits are made separately for 
each intensity bin. By using spatially high-pass filtered measurements of the LOS component of the magnetic field, the flow properties 
are determined separately for magnetic spines (relatively strong and vertical field) and inter- spines (weaker and more horizontal field). 
Results. The dark convective downflows discovered recently in the 538.03 nm line are evident also in the 630.15 nm line, and have 
similar strength. This convective signature is the same in spines and inter- spines. However, the strong radial (Evershed) outflows are 
found only in the inter- spines. 

Conclusions. At the spatial resolution of the present SST/CRISP data, the small-scale intensity pattern seen in continuum images 
is strongly related to a convective up/down flow pattern that exists everywhere in the penumbra. Earlier failures to detect the dark 
convective downflows in the interior penumbra can be explained by inadequate spatial resolution in the observed data. 

Key words. Sunspots - Convection - Magnetic fields - Magnetohydrodynamics (MHD) 



1. Introduction 

The origin of sunspot penumbral filam entary structure and 
strong horizontal flow (lEvershedl 1 19091) . their complex mag- 
netic field topology, and the mechanism that transports en- 
ergy to the visible surface represent some of the most long- 
standing puzzles in astrophysics. An early proposal for explain- 
ing this filamentar y structure was in terms of convection rolls 
(iDanielsonll 1 96 1 bh . based on sun spot photographs r ecorded dur- 
ing the 1959 Stratoscope flights (iDanielsonll 1 96 1 ah . Initial sup- 
port for this interpretation came from spectroscopic observa- 
tions at l ow spatial resolution of a sunspot at 16° heliocentric 
distance (iBeckers & Schrote^ 1 19691) . demonstrating a correla- 
tion of 29-39% between the continuum intensity and LOS ve- 
locity measured in the line wings of the Fei line at 557.6 nm. 
Later observations did not give consistent suppor t for strong 
such c orrelations (e - g-- Iwiehr & Stellmache3 Il989l: iLites et aP 
Il990l) . iJohannessonI (1 19931) found a generally weak or non- 
existent correlation between the continuum intensity and the 
flow field, with only a tendency for bright structures to cor- 
respond to local blue s hifts o n both sides of the spot, while 
ISanchez Almeida et al.l (1 19931) claimed clear evidence of such 
a correlation, based on 21 slit Fei 630.15 nm and 630.25 nm 
spectra recorded with an 0'.'5 spectrograph slit. More recent ob- 
servations with the Swedish 1-m Solar Telescope (SST) at much 
higher spati al resolution of part of a sun spot at 35° heliocen- 
tric distance (ISanchez Almeida et aD[2007l) demonstrated a clear 
correlation between intensity and vertical velocity in the sense 



expected for convection. We note, however, that none of the 
above observations demonstrated evidence that the dark penum- 
bral filaments actually have absolute downward flows, as re- 
quired for these flows to be identified as convective. We finally 
note that indirect evidence for penumbral convection was found 
from analysis of proper motions of small-scale intensity struc- 
tures, obtained by employing local cro ss correlation techniq ues 
to a time series of SST G-band images (iMarquez et al .1120061) . A 
relevant question is to what extent velocities inferred from inten- 
sity patterns reflect true flows. 

In parallell, attempts were made to explain the Evershed 
flow, and to characterize its relation to bright and dark fila- 
mentary fine structure. A highly infl uential explanation was i n 
terms of siphon flows in flux tubes (iMeyer & Schmidtl 1 19681) . 
This apparently triggered the development of models and nu- 
meric al simulations of flows in flux tubes for a period of over 30 
vears (iThomas & Montesinoslll993l:ISolanki & Montavon|[T993l: 



ISchlichenmaier et al.lll998l: ISchlichenmaieill2002l) . Support for 
such models came from observations at about O'.'T spatial reso- 
lution showing bright upflows in th e inner penumbra and dark 
downflows in the outer penumbra ( Schlichenmaier & Schmidj 
119991: ISchmidt & Schlichenmaied l2000l) . These authors also 
found an overall correlation between brightness and vertical 
velocity that can be interpreted as an indication of convec- 
tion. However, throughout the inner two thirds of the penum- 
bra observed, both the locally bright and dark component of the 
penumbra on the average showed only upflows. 
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An important step toward our understanding of the 
fine structure of penumbral magnetic fields came with the 
discovery of broa dband circular polarization in sunspots 
(lining et alJ Il974allbb . implying the simultaneous existence of 
strong line-of-sight (LOS) gradients in the magn etic field 
and flow velocity (ISanchez Almeida & LitesI Il992h . The in- 
ferred strong gradients appeared to imply extremely strong 
currents and curvature f orces in the vi si ble layers of penum- 
bra^_^anchez Almeida & Lites| Il992t ISolanki et al] 119931: 
ISolanki & Montavonii 19931) (note however, that strong gradients 
in the magnetic field do not necessarily imply strong forces: po- 
tential fields above nearly field-free gaps giv e rise to strong LOS 
gradi ents above the gaps, but are force-free : ISpruit & Scharme^ 
I2OO6I) . This (appa rent) problem was addressed by the uncombed 
penumbra model (ISolanki & Montavonll 19931) . based on the as- 
sumed existence of slender nearly horizontal flux tubes with con- 
stant magnetic field within the flux tube (no volume currents). 
The flux tubes were assumed to be elevated above the penumbral 
photosphere and carrying the Evershed flow, and to be embedded 
in a more vertical background magnetic field. V arious calcula- 
tions andjnyersions^ & Mont avonll 1 9931: [M artinez Fillet 
200Q[ ISchHche nm"^r & Collado si 120021: iBellot Rubio et al. 
20031 12004: Borrero et al.. 2007: iTritschler et al.l l2007h based 
on this and similar simplified flux tube models gave consis- 
tency with the spectropolarimetric data, at that time rarely 
reaching better than V spatial resolution. An alternative ex- 
planation of obse r ved penumb ral Stokes spectra was given by 
I Sanchez Almeidal (l2005l l2006l) . He postulated the existence of 
micro- structured magnetic field and velocity structures (repre- 
sented as two opticall y thin magnetic compon ents) and used a 
special inversion code (ISanchez Almeidairi997l) to test this con- 
jecture. Fitting his observed spectra required the existence of si- 
multaneous upflows and downflows throughout the penumbra, 
which can be interpreted as evidence of convection. However, 
these inversions were constrained by the temperature being equal 
for the two magnetic components, which excludes any convec- 
tive energy flux. A surprising result of the inversions is the exis- 
tence of opposite polarity field everywhere in the penumbra. We 
note that the inversions constrain the magnetic fields and flows 
in the two magnetic components to be strictly parallell, such that 
any mixture of upflows and downflows must be associated with 
mixed polarity field, and vice versa. 

It was soon evident that a flow in a slender flux tube 
cannot carry the energy needed to compensate the energy 
losses throughout the en t ire ra dial extent of a penumbra 
(ISchlichenmaier & Solankil l2003l) . This a nd other problems 
of flux tube models were high-lighted by ISpruit & Scharmeil 
(I2OO6I) . They concluded that the origin of the filamentary struc- 
ture and the large fluctuations in field strength and magnetic field 
inclination across penumbral filaments is not embedded flux 
tubes but must be convection in (nearly) field-free gaps below 
the visible surfaceQ, resolving the heat flux p roblem. Simple po- 
tentia l field and magnetohydrostatic models ('Scharmer & Spruij 
[2OO6I) of such "gaps" demonstrated that large fluctuations in 
magnetic field strength and inclination are unavoidable conse- 




^ We make two remarks. First, as in the quiet Sun, the penumbral 
photosphere is convectively stable (lRempelll2011l) . and the flows seen 
in the line forming layers represent the overshoot of this convection. 
Gaps are thus expected primarily in the invisible layers just below 
the photosphere, and to what extent they are field-free can probably 
only be investigated by means of simulations. Second, the gaps are 
not expected to be strictly field-free, but rather to have sufficiently re- 
duced field strengths to allow the overturning co nvective flows to occur 
(iNordlund & Scharmeiil20TQl : ISpruit et al.ll20TQl) . 



Fig. 1: Overview of the entire FOV observed. The highlighted 
contiguous area close to the boundary of the FOV corresponds 
to the area used as quiet Sun reference, the highlighted area in 
the penumbra corresponds to what is referred to as the "interior 
penumbra". The top panel corresponds to the 538 nm continuum, 
the right bottom to the COG velocity measured in the same line. 
Tick marks are at intervals. The arrow points in the direction 
of Sun center. The dark rectangles outline sub-fields shown in 
Figs.Oandlll 



quences of such gaps, but without giving rise to any magnetic 
forces, and that the magnetic field should be much more hori- 
zontal above than outside the gaps. The model also explained the 
penumbral dark cores ( Scharmer et al. 2002) as due to a strongly 
warped t= 1 surface caused by the strong azimuthal fluctuations 
in gas pressure (induced by the strong fluctuations in magnetic 
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Fig. 2: The panel shows LOS component of the magnetic field 
over the same FOV as in Fig. [T] Tick marks are at V intervals. 
The arrow points in the direction of Sun center. 



field strength) combined with an overall drop in temperature 
with heighQ. While the convective gap model did not at that time 
explain the origin of the Evershed flow, it was evident that this 
flow can only exist in the convective gaps or immediately above 
them , where the field is nearly horizontal (I Spruit & Scharmerl 
120061) . 

Theoretical support for the convective origin of penumbral 
fine structure came with the first MHD sim ulation of penum- 
bral fine structure (iHeinemann et al.l l2007l) . While limited to 
a thin slice of a quite small modelled sunspot, these simula- 
tions demonstrated convection with strongly reduced (but far 
from zero) field strength, a systematic (but weak) radial out- 
flow in gaps, dark cores above short penumbral filaments, an in- 
ward migration of the penumbral structure toward the penumbra 
(as observed in movies), and a moat flow outside the sunspot. 
Based on these simulations, the Evershed flow was explained 
as being id entical to the horiz ontal component of penumbral 
convection (IScharmer et al.|[20Q 8b). Sim ulations of miich larger 
sunspots and w ith higher grid resolution (iRempel et al.ll2009allbl : 
lRempell[2QTTh demonstrated a higher degree of realism, with 
much longer filaments and stronger Evershed flows. Detailed 
analysis of the simulation data has clarified some of the prop- 
erties and drivi ng mechanisms of the penumbral convection and 
Evershed flow (lRempelll201ll) . Remarkably, the length scales, 
mass fluxes and RMS velocities of penumbral convective flows 
are quite sim ilar to those of field-free convection in the quiet Sun 
(iRempel et al . 2009a). 

Strong evidence for penumbral convection came from obser- 
vations of filame ntary dynamics i n sunsp ots located well away 
from disk center. Ilchimoto et al.l (l2007bh found "twisting mo- 
tions" in filaments from space-time plots along lines crossing 

^ This is very similar to what is seen in faculae and bright points 
where locations of strong magnetic field correspond to a locally de- 
pressed T = 1 surface, leading to enhanced intens ity due to the higher 
temperatures in the deeper layers jScharmedl2009l) . 



filaments in the inner penumbra. Further evidence for convec- 
tion (interpreted as convective rolls) in penumbral filaments 
in SST observa tions of a sunspot at 40° heliocentric distance 
was reported by lZakharov et al.l (l2008l) . More recently, evidence 
for penumbral convectio n was observed in a sunspot located at 
58° heliocentric distance (ISpruit et al.ll2010l) . These authors also 
pointed out that the tilted "striations", seen in the filaments are 
similar to those observed in fa culae close to the lim b and re- 
produced in MHD simulations (ICarlsson et al.ll2004l) . and most 
likely outline magnetic field lines. 

Whereas the convective origin of the penumbral structure 
and Evershed flow thus seems quite clear from a theoretical 
perspective, and also was receiving strong support from obser- 
vations, the accumulating evidence for penumbral convection 
was not readily accepted. The main objection is the absence of 
direct observ ational evidence for convective dow nflows inside 
the penumbra. iFranz & Schlichenmaie j (l2009l) and lFranj (1201 ll) 
analyzed spectropolarimetric data from Hinode, recorded from 
sunspots close to disk center and found no indication s of con- 
vection in penumbral filaments. However jFranzl (1201 ll) actually 
measured a correlation coefficient for small-scale (high-pass fil- 
tered) intensity and velocity fluctuations of up to -53% in the 
inner parts of one penumbra in the Fei 630.15 nm line (here- 
after referred to as "6301 line"), but considered this degree of 
correlation as being too small to be significant. 

While clear evidence for convection thus appears to be miss- 
ing (or, is claimed to be missing) in spectropolarimetric Hinode 
data, the Swedish 1-m Solar Telescope (SST) pro vides clear 
such evidence. We recently (Sch armer et al.l l2Qll1) made ob- 
servations in the Ci 538.03 nm line (hereafter referred to as 
the "5380 line"), formed very clo se to the visible photosphere 
(ISchlichenmaier & Schmidtl 1 19991) . where (overshooting) con- 
vective flows are expected to be significantly stronger than at 
the heights of formation of the 6301 line. The spatial resolution 
of this SST data is more than 2.3 times higher than that of the 
spectropolarimetric Hinode data, representing a major improve- 
ment in our ability to resolve penumbral filamentary structure 
and flows. By making azimuthal fits to separately determine the 
variation of radial (horizontal) and vertical velocities with inten- 
sity, clear evidence of a correlation between intensity and verti- 
cal velocity was found in the sense expected for convection. In 
addition, the darkest penumbral features in the interior penum- 
brc0 show downflows up to about 1 km s~^ while the brightes t 
features show upflows of several kms"^ (IScharmer et al.l l20Tl1) . 
Additional support for the existence of convective downflows 
was recently found fr om analysis of other SST data in the 5380 
line (iJoshi et al.ll201ll) . Usin g simulated SST obser v ations based 
on num erical simulations of IRempel et al.l (l2009b). Bh arti et aD 
(1201 ll) recently confirmed that convective downflows should 
be observable with the 5380 line at the spatial resolution of 
the SST, while such downflows would be more difficult to ob- 
serve in the Fei 709.1 nm line, explaining in part the failure of 
iBellot Rubio et al. l (l20Tol) in detecting such flows in this line. 
The other diflSculty in detecting vertical flows from their data is 
the contamination of the LOS velocity from the stro ng horizontal 
Evers hed flows, due to the 5? 4 heliocentric distance (iBharti et al.l 

Honl). 

In the presen t paper we extend the analysis of the observa- 
tions reported by IScharmer et al.l (l201ll) to include Doppler and 
magnetic field measurements made in the 6301 line, recorded 



^ Here and in the following we use the term "interior penumbra" to 
mean the part of the penumbra that is well away from both the inner 
(umbral) and outer (quiet Sun) boundaries. 
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nearly simultaneously with the recently published 5380 data. We 
demonstrate that the convective signatures found in the 5380 line 
also are obvious in the wings of the 6301 line, although with 
reduced velocity amplitudes, and that the dark penumbral fea- 
tures in the interior penumbra on the average are associated with 
(convective) downflows also when observed in the 6301 line. We 
also find that the convec tive sig natures are the same in magnetic 
spines and inter- spines (iLites et al. 1993), but the strong radial 
outflows are seen only in the inter-spines. In Sect. 2 of this pa- 
per, we describe the observations, data processing and methods 
of data analysis. Correlations between continuum intensities and 
radial and vertical velocities obtained from azimuthal fits of LOS 
velocities measured in the 5380 and 6301 lines are compared and 
discussed in Sect. 3. Finally, in Sect. 4 we discuss the results and 
draw conclusions with respect to models. 

2. Observations and data processing 



The data were recorded with the SST (IScharmer et al.ll2003ah 
using the CRisp I maging SpectroPolarimeter (CRISP: IScharmeij 
l2006t IScharmer et a l. 2008a). The SST is an un-obscured 1- 
meter aperture telescope th at presently employs lo w-order (37 
electrode) adaptive optics (IScharmer et al.l l2003bh to achieve 
nearly diff'raction limited performance in excellent seeing at 
wavelengths as short as 390 nm. Of importance in explaining 
the excellent data quality is the use of multiple short (17 msec) 
exposures that are p ost-processed with an e laborate image recon- 
struction technique (Ivan Noort et al.ll2005l) to compensate resid- 
ual low-order aberrations left by the AO system. CRISP is a dual 
Fabry-Perot filter system with the low-resolution etalon having 
reduced reflectivity in order to mitigate the efl'ects of cavity er- 
rors (variations in separation between the reflecting surfaces of 
the etalons). The telecentric mounting of the etalons gives wave- 
length shifts of the transmission peak over the FOV, but the low 
reflectivity of the low-resolution etalon gives a passband that is 
wide enough to accommodate the relative wavelength shifts be- 
tween the two etalons, such that the wavelength transmission 
profile shows only small variations in its shape over the FOV. 
The spectral resolution of CRISP is modest, about 4.4 pm at 
538 nm and 6 pm at 630 nm, in order to allow fast tuning of 
spectral lines at high signal-to-noise and without spectral un- 
dersampling. To support image reconstruction and co-alignment 
between images recorded at diff'erent wavelengths, the setup in- 
cludes one CCD camera that records broad-band images through 
the CRISP pre-filter used to select a particular spectral line. To 
enable polarimetric measurements with low levels of seeing- 
induced cross-talk from Stokes / to Q, U, and V there are two 
cameras at the CRISP final focal plane, with the light divided be- 
tween the two cameras by a polarizing beam splitter. Polarization 
modulation is made with two tunable liquid crystals (LC's) with 
their fast axes at 0° and 45° with respect to the vertical axis on 
the optical table. The LC's have sufficient stroke to allow full 
Stokes measurements over the entire wavelength range of CRISP 
(from about 510 nm to 860 nm). 

The present data were recorded on 23 May 2010 at 14: 1 1 UT 
from a reasonably regular sunspot at 15° heliocentric distance. 
The data consists of three nearly simultaneous (the time differ- 
ence between the 5380 and 6301 line scans is 13 s) spectral scans 
of the 5380 line, the 6301 line, and the Fei at 630.25 nm (here- 
after the "6302 line"). Due to the telluric blend in the red wing 
of the 6302 line preventing measurements of the bisector wave- 
length shifts close to the continuum, we will not discuss the 6302 
data in the present paper. The images were recorded at an image 
scale of (y.'059 per pixel. After image processing, co-alignment 



of the images recorded at the three diff'erent wavelengths and re- 
moving the apodized part of the FOV near the edges, the science 
data consisted of 881x881 pixels covering 52''x52'', shown in 
Figs. [T] and [21 The 5380 data consist of a total of 960 images 
sampled with CRISP at 23 wavelengths in typical steps of 2.8 pm 
plus 480 broad-band images. The 6301 (and 6302) data consist 
of a total of 600 images sampled at 15 wavelengths in steps of 
4.4 pm and 320 broadband images. 

Figure [T] shows the 538 nm continuum image (top) and the 
5380 center-of gravity (COG) line-of-sight (LOS) velocity of the 
entire field-of-view (FOV) observed. Fig. [2] shows the LOS mag- 
netic field ob tained with the COG method (iRees & Semelll 19791: 
ICauzzi et al. The high-lighted region around most of the 

outer boundary of the FOV defines the area used as quiet Sun 
reference (although about a few per cent of this area obviously 
is covered with strong field), the high-lighted area in the sunspot 
outlines what we refer to as the "interior penumbra", which is 
the main target of our analysis. Figures [3] and |4] show all mea- 
sured quantities within sub-fields (marked with dark rectangles 
in Fig.[T]) close to the limb-side and disk center-side parts of the 
penumbra. 

We note that Fig.[2]shows very little evidence of the opposite 
polarity fields observed well inside the penumbra in a sunspot 
located at small heliocentric distances (llchimoto et al ] 12007^ 
lFranzll201ll) . However, Fig. |4] (upper row, mid panel) shows sev- 
eral small-scale patches where the LOS magnetic field is very 
weak at the center-side penumbra. The absence of opposite po- 
larity patches in our LOS magnetic map may be related to the 
associated strong downward ffows and the corresponding abnor- 
mal Stokes V profiles that are not well captured with our simple 
(COG) method for estimating the LOS magnetic field. However, 
inspection of the Stokes V images at -h26.3 pm and -h30.7 pm 
confirms that opposite polarity patches are found primarily in 
the outer penumbra and that only a few such patches are seen 
well inside the penumbra. 

2.1. MOMFBD image reconstruction and alignment 

The images were corrected for gain and offset, using flats at a 
fixed continuum wavelength. The gain and dark corrected im- 
ages from the 5380 and 6301/6302 lines and their correspond- 
ing broadband images were processed separately as two data 
sets with the Multi-Obje ct Multi-frame B lind Deconvolution 
(MOMFBD) method ( van Noort et al.ll2005h . In the MOMFBD 
image reconstruction process we conservatively compensated 
for only 36 Karhunen-Loeve aberrations. This leaves the small- 
scale wavefront errors uncompensated for at scales roughly 15 
cm or smaller. With sufficiently good seeing, as in our case, the 
effect of such truncation is to leave residual straylight in the 
wings of the point- spread-function (PSF), whereas the core is 
well corrected, but wit h reduced peak strength (red uced Strehl) 
(Sc harmer et"aDl2010l: iLofdahl & Scharmed [20lTh . We expect 
the restored images to have a spatial resolution close to the 
diffraction limit of Off 14 at 538 nm and a' 16 at 630 nm. 

After image reconstruction, the recorded Stokes images were 
demodulated for telescope polarization using the SST polariza- 
tion model developed by Selbing ( 2005) Finally, the 6301 and 
and 5380 data were co-aligned by destretching the 6301 contin- 
uum image onto the 5380 continuum image over 16x16 sub- 
fields. The so-obtained matrix of geometrical distortions was 
then applied to all measured 6301 quantities. This process com- 
pensates to some extent also for any small effects of large scale 
flows and evolution of the structure during the 13 sec time 
elapsed between the recording of the 5380 and 6301 data. 
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Fig. 3: The Figure shows a 5''x8'' sub-field from the right-hand side (mostly limb-side) penumbra, outlined in Fig.[T] The top row 
shows the 538 nm continuum, the 5380 COG velocity within the spine mask, the LOS magnetic field, the 5380 COG velocity within 
the inter- spine mask, and the 6301 continuum intensity. The bottom row shows the 5380 COG velocity, the 5380 line core velocity, 
the 6301 COG velocity, 6301 70% bisector velocity and the 6301 line core velocity. The velocity maps have been scaled individually 
to enhance fine structure. Maps shown represent original quantities and have not been processed with unsharp masking. Tick marks 
and small white dots are at V intervals. 



2.2. Cavity error compensation and wavelength calibration 

In principle, it would be sufl&cient to measure the wavelength 
shifts from the cavity errors and to use these to correct the mea- 
sured LOS velocities. However, stray light compensation of the 
science data also redistributes spatially the wavelength shifts 
from any cavity errors. Therefore, we first shifted the line pro- 
files to a common wavelength scale by compensating for the cav- 
ity errors, and then applied the stray light compensation to the 
data. 

The cavity error wavelength shifts were measured from the 
averaged flat-field images by fitting second-order polynomials 
to the core of the Stokes I profile for each line. The wavelength 
shifts were determined from the coefficients of the fits. The 
broad 5380 line is well sampled (at 2.8 pm wavelength steps) 
and we used linear interpolation to shift the line profiles to a 
common wavelength scale and at the same time applying a first- 
order correction to compensate for the transmission profile of the 
CRISP pre-filters. The 6301 line has coarser wavelength sam- 
pling (4.4 pm) and we used cubic splines to resample the line 
profiles to 1 . 1 pm steps before shifting the line profiles with lin- 
ear interpolation according to the measured cavity errors. After 
compensation for cavity errors, we calculated average 5380 
and 6301 line profiles from the high-lighted region surrounding 
most of the sunspot, shown in Figs. [T] and O We calibrated the 
wavelength scales such the line core velocities measured from 



the spatially av eraged 5380 and 6301 line profiles correspond 
to -917 ms"^ (IScharmer et al.1 l201lh resp. -200 ms"^ (neg- 
ative velocities correspond to blue- shifted profiles), in agree- 
ment with the convective blue-shifts obtained from simula- 
tions (Ide la Cruz Rodriguez et al.l [201 ll) . The 6301 convective 
blue- shift adopted is somewhat larger than the shift obtained 
from convection simulations (-188 ms"^) at disk center for 
the line core, when the averaged qui et Sun profile is convolved 
with the CRISP transmission profile (Ide la Cruz Rodriguez et al.l 
[mil). As a check of this calibration, we calculated the aver- 
age quiet Sun LOS velocity from the spatially resolved spec- 
tra and obtained -130 ms"^ and -220 ms"^ for the COG 
and line core 5380 velocities, and 150 ms"\ 120 ms"^ and 
-100 ms"^ for the 6301 COG, 70% bisector and line core ve- 
locities. Averaging over the entire penumbra, the corresponding 
values are -200 ms"^ and -290 ms"^ for the 5380 line, and 
50ms"^ 30ms"^ and -130 ms"^ for the 6301 line. The remain- 
ing convective blue-shifts in the 5380 line are similar to those ex- 
pected from convection simulations (see SOM, Table S2). The 
values for the 6301 line are reasonably consistent and indicate 
that the systematic errors are not larger than ±150 m s"^ for this 
line. The weak blend in the red wing of the 6301 line most likely 
has a small influence on the velocities measured from the wings 
(the COG and 70% bisector measurements) of the 6301 line, but 
we made no attempt to compensate for that. 
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Fig. 4: The panel shows a 5''x8'' sub-field from the left-hand side (mostly disk center-side) penumbra outlined in Fig. [T] The 
arrangement of the panels is the same as shown in Fig.[3l 



The sunspot is supposedly surrounded by a r adially outward 
moat flow (e.g.. I Vargas Dominguez et al]|2QQ8l) of up to about 
0.3-0.5 km s"^ that will have an influence on our wavelength cal- 
ibration. No obvious evidence of such a flow can be seen in our 
Doppler maps, even with heavy spatial smearing of the individ- 
ual Doppler maps. Since the region used for wavelength calibra- 
tion surrounds most of the sunspot and because of the relatively 
small heliocentric distance, we expect the systematic errors in 
the wavelength calibration from the moat flow to be small. 

2.3. Stray light calibration and compensation 

Recently, the RMS contrast of solar granulation was measured 
from Hinode data, combined with a detailed characterization 
of the str aylight properties o f the t elescope and instrumenta- 
tion used. I Wedemey er-B ohml (l2008l) determined the straylight 
PSF from images recorded during both a Mercury transit and a 
partial solar eclipse. After deconvolution of the granulation im- 
ages for the measured stray-light PSF, the corrected granulation 
contrasts are close to those obtained from two independently 
developed MHD si mulation codes (the CO^BOLD and Stein- 
Nord l und codes; |W edemev er-Bohm & Rouppe van der VoortI 
l2009l) . IPanilovic et al . ( 2008) degraded synthetic continuum im- 
ages at 630 nm calculated from simulations with the MURAM 
code with the PSF of Hinode. Taking into account only the PSF 
of a perfect unobscured 50 cm telescope, the RMS contrast is 
reduced from 14.4% to 10.9%. Including the efl'ects of the large 
central obscuration and spider plus the CCD, reduces the RMS 
contrast to 8.5%. This is still higher than the observed value of 
7%. Danilovic et al. interpret this difl'erence as the result primar- 



ily of small focu s errors and/or other residual aberrations, bu t 
the work of W edemey er-B ohm & Rouppe van der VoortI (l2009l) 
as well as the Gaussian fits of iMathew et a" ■I (I2009h clearly in- 
dicate significant contributions also from the far wings of the 
PSF. We conclude that there are no longer scientific grounds 
for questioning the RMS granulation contrast values obtained 
from numerical 3D simulations, and that the measured granu- 
lation contrast in highly resolved granulation images recorded 
with the SST and other solar telescopes therefore can be used to 
constrain the amount of straylight. We also note that the large 
central obscuration, the wide spider and aberrations and/or other 
imperfections combine to reduce the observed RMS contrast of 
granulation to about 64% of that theoretically possible at wave- 
lengths around 630 nm for a telescope with the diameter (but 
without the spider and central obscuration) of Hinode. 

For the SST, the observed RMS contrast is around 6.5-7% 
at 630 nm in the best "raw" images, but since the spatial res- 
olution is higher than for Hinode, we could in principle reach 
14% contrast with a perfect SST operating outside the Earth's at- 
mosphere. MOMFBD image reconstruction increases the RMS 
contrast somewhat to about 7.5% at 630 nm and 9% at 538 nm 
(IScharmer et al.l 120101) . Clearly, the SST data must contain on 
the order of 50% straylight, operating at scales correspond- 
ing to granulation or larger. At the same time, "conventional" 
straylight, often characterized by a PSF having very wide (e.g., 
Lorentzian) wings, must be low in this and other SST data: the 
minimum umbral intensity in the sunspot presently discussed is 
15% at 538 nm and 18% at 630 nm. Measurements of straylight 
made during the 2004 Venus transit suggests on the order of 4% 
at 390 nm and 9% at 700 nm (lKiselmanll2008l) . Measurements of 
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stray light from the primary focal plane just outside the SST fo- 
cal plane and including all relay optics as well as CRISP and the 
CCD itself sugg ests levels of conventional" straylight of at most 
a few percent (iLofdahl & Scharme3 1201 ll) . The origin of the 
dominant source of straylight is largely unknown but is now sus- 
pected to be mostly from aberrations, including the high-order 
aberrations truncated in MOMFBD processing, high-altitude 
seeing and/or effects from the rela t ively long integration time s 
(17 ms) used (IScharmer et aDl201Ql:lL5fdahl & Scharme3l201lh . 
This conclusion receives at least partial support from a recent 
comparison between speckle and MOMFBD image restorations 
of spectropolarimetric data recorded at the German VTT. This 
comparison demonstrated that speckle restorations and MFBD 
processing (identical to MOMFBD processing, but with a single 
wavelength channel) gives significantly higher granulation con- 
trast than with MOMFBD (iPuschmann & Beck 201 1). 

An important component of our data processing is 
therefore straylight compensatio n. In th e Supplementary 
Online Material (hereafter SOM) of lScharnie r et al. (201 1) (see, 



be very wide. In particular, a PSF having extended Lorentzian 
wings would lead to observed umbra intensities on the order of 
50%, given that the level of straylight must be on the order of 
50% to explain the observed granulation contrast. If we assume 
that the true minimum umbra intensity is not less than half of 
the measured values, then up to 7-9% of the straylight could be 
from a Lorentzian PSF. A similar problem with the umbra inten- 
sity arises even with a Gaussian PSF, if it is too wide. Based on 
extensive tests, we conclude that the majority of the straylight 
must be contained inside a PSF with a FWHM that is less than 
about 1" else the restored um bra intensity is close to zero or 
negative (IScharmer et al.|[20Tl]) . 

Our best agreement with granulation data, taking into ac- 
count the constraints from the umbra intensity, is a Gaussian 
PSF with a FWHM of about only V'.2 and a straylight con- 
tribution of a=58% at 538 nm0 We have adopted these val- 
ues for the 5380 line. We use the same FWHM for the 6301 
line, but with lower straylight contribution, Qf=50% at 630 nm, 
since our present observations as well as those discussed by 



http : //www . sciencemag . org/content/3 3 3/6040/3 16/ supri^^ a(etiever-Bohm & Rouppe van der VoortI (l2009l) suggest that 



we discuss this in detail for the 5380 line, here we summarize 
the analysis made. W e used snapshots from fi eld-free simula- 
tions obtained with the lStein & Nordlundl(ll998l) code, degraded 
with the PSF of a perfect 1-m solar telescope. This constitutes 
our "true" theoretical data. We assumed that the actual PSF of 
the MOMFBD restored science data corresponds to a diff'raction 
limited core with reduced peak strength (reduced Strehl) and 
enhanced wings due to straylight as follows 



lo = (l-a)It + aIt^S(W), 



(1) 



where lo and It are the observed and "true" intensities at any 
wavelength and polarization state, a the straylight fraction, "*" 
denotes convolution and S is the straylight PSF, having a full 
width at half maximum (FWHM) W. This corresponds to a total 
PSF P of 

P = (l-a)S-^aS(W), (2) 

where S is Dirac delta function. Note that this PSF corresponds 
to assuming that the core of the PSF is perfectly compensated for 
by MOMFBD processing; this is consistent with our approach 
of first degrading the synthetic simulation data with the diff'rac- 
tion limited PSF of the SST. When using Eq. (2) to deconvolve 
the data we thus compensate only for straylight. We used this 
PSF to both degrade the simulated data and to deconvolve the 
observed data. As basis for our comparison, we calculated both 
synthetic continuum images and line profiles corresponding to a 
heliocentric distance of 15°. We then calculated 5380 line core 
LOS velocities from both the synthetic and the observed data. 
The calculated quantities compared were the RMS continuum 
intensity, the 5380 line core RMS velocity and the spatially av- 
eraged line core velocit}0. This comparison does not give much 
guidance about the shape of the straylight PSF. However, the low 
umbra intensity observed implies that the straylight PSF cannot 



^ The Ci 5380 line is very temperature sensitive and weakens con- 
siderably in dark granular and penumbral lanes. This makes it very 
sensitive to straylight in dark (cool) lanes, explaining the large con- 
vective blueshift (-917 ms~^ at 15° heliocentric distance). This also 
makes the spatially averaged 5380 velocity from spatially resolved line 
profiles a sensitive indicator of straylight. Indeed, even the straylight 
compensated SST data leaves a residual convective blueshift of about 
-240 m s~^ when averaging LOS velocities from the spatially resolved 
line profiles for quiet Sun (this residual is consistent with what is ob- 
tained from simulations). A similar blue- shift is obtained when averag- 
ing the 5380 line core velocities over the penumbra. 



the SST straylight increases at shorter wavelengths. This leads to 
a quiet Sun RMS granulation contrast of about 1 3.3% at 630nm, 
whic h is somewhat lower th an from simulations (iDanilovic et al.l 
2008 HScharmer et al.ll2010l) . but the "quiet" region we observed 
does contain some small fraction of strong field that is likely 
to reduce the overall contrast. Admittedly, this determination of 
the straylight PSF is not precise, but the main eff'ect of the stray- 
light is present already when its PSF includes about one granule. 
Compensation with a wider straylight PSF will require some- 
what smaller values of a, but the overall results do not change 
qualitatively. Whereas the accuracy of our measured intensities 
and LOS velocities thus is limited by the uncertainties about the 
precise shape of the straylight PSF, we are confident that the con- 
clusions drawn about the clear existence of correlations between 
measured continuum intensities and LOS velocities, as well as 
about the existence of dark convective downflows in the interior 
penumbra, are robust. 

2.4. LOS velocities and magnetic field measurements 

When measuring spectral line properties at low temperatures 
(low continuum intensities) there is always a concern about the 
possible influence of blends from (mostly) molecular lines. The 
5380 line weakens considerably at reduced temperatures and 
must disappear in the umbra, strongly aggravating any such in- 
fluence of blends. In SOM, we established that the influence of 
blends is small near the core of the 5380 line at continuum in- 
tensities corresponding to granulation and penumbral filaments. 
However, there is evidence of weak blends both in the far blue 
and red wings (see Fig. S5 in SOM) in the penumbral dark struc- 
tures and in the dark intergranular lanes. These blends are obvi- 
ous in the umbra (where also a strong molecular blend, shifted by 
about 3-3.5 pm to the blue of the C i line, is present - see SO1V0). 
To mitigate the influence of the blends in the wings, we initially 
choose to estimate LOS velocities from the line core only, by 
fitting a second -order polynomial to 5 wavelengths closest to the 
line minimum (IScharmer et al.ll2QTl1) . In the present work, we 
estimate also 5380 COG velocities, expected to be dominated 



^ A narrow straylight PSF is consistent with main contributions from 
high-order aberrations, such as fr om seeing or possibly from the adap- 
tive mirror ( Scharme r et a l. 201Ql : lLofdahl&Scharmer 201 L) . 

^ This blend has now been identified to be from Mg H by Uitenbroek, 
Dumont and Tritschler (in prep.). 
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Fig. 5: Examples of measured quantities with superimposed contours outlining the six radial zones in the penumbra, numbered 1 
(adjacent to the umbra) to 6 (the outermost penumbra). Radial zones 2-4 constitute what is referred to as the "interior penumbra". 
In the top row are shown (left to right) the 538 nm continuum intensity, the COG velocity in the same line and the LOS component 
of the magnetic field, obtained with the COG method. In the lower row are shown the same quantities but high-pass spatially filtered 
to remove large-scale azimuthal and radial variations. Tick marks are at intervals. The arrow points in the direction of Sun center. 



by flows even closer to the photosphere than from the line core. 
In order to remove influence from the above mentioned blue/red 
blends, we crop the 5380 line profiles by excluding the first 6 and 
last 2 wavelengths shown in Fig. S3 in SOM. This leads to un- 
derestimates of the strengths of the strongest blue- and redshifted 
COG velocities, as is the case when the COG method is ap plied 
to LOS magnetic field measurements (ICauzzi et al.lll993l) . but 
these data nevertheless add support to the earlier line core mea- 
surements and also provide interesting information about flows 
in the deepest observable layers of the penumbra. 



We measured the 6301 line core shifts from the 3 wave- 
lengths nearest to the line minimum, the COG velocity from 
the entire line profile and (using the spline interpolated line 
profiles) the 70% bisector velocity, where the 70% intensity 
level is relative to the line minimum intensity at each pixel. 
We estimated the LOS magnetic field ^ lqs by applying the 
COG method to Stokes 7 -h V and 7 - V (iRees & Semell 1 19791: 
ICauzzi etal ] ll993l:IOrozco Suarez et al.ll2010l) and from the dif- 
ference in the two COG wavelengths SA we obtained ^los = 
SA (9.334 10~^^ A^gi)'^, where wavelengths are given in units of 
Angstrom, and gi is the Lande splitting (=1.667 for the 6301 
line). 



2.5. High-pass spatial filtering and masks 

In order to investigate relations between small-scale fluctuations 
in intensity, LOS velocity and magnetic field, we make repeated 
use of high-pass filtered (unsharp masked) quantities. In particu- 
lar, such filtering strongly reduces radial variations of the inten- 
sity in the penumbra, large-scale radial and azimuthal variations 
in the LOS component velocities from the horizontal (Evershed) 
flow and the horizontal component of the magnetic field. For 
example, the high-pass filtered continuum intensity Sic is calcu- 
lated as 

SIc = Ic-G(W2)^Ic (3) 

where G is a Gaussian profile having a FWHM width W2 of 
20 pixels, or approximately V!2. The same high-pass filter is 
used with the measured LOS velocities and LOS magnetic field. 
The FWHM used here is coincidentally the same as used for the 
straylight compensation, but the precise value of the FWHM is 
not critical for the analysis. It might be argued that a filter kernel 
with a larger FWHM would be more appropriate, but we prefer a 
narrow filter in order to reduce the spatial mixing of information 
across the borders to the umbra and the surrounding quiet Sun. 

An obvious eff'ect of the spatial filtering is that the abso- 
lute zero point of any filtered quantity is lost and that all fil- 
tered quantities are referenced only to a local mean, defined by 
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the filter kernel. In the following, we therefore analyze both the 
filtered velocities to show their relation to local intensity fluc- 
tuations, and the absolute (unfiltered) velocities to demonstrate 
that the darkest penumbral structures are indeed associated with 
downflows. 

In Fig. [5] are shown in the top row the 538 nm continuum 
image, the 5380 COG velocity and the LOS component of the 
magnetic field. In the bottom row are shown the same quantities 
after applying the previously defined high-pass filter. Clearly, 
the high-pass filter strongly reduces both systematic radial vari- 
ations as well as diff'erences between the limb and center sides 
of the penumbra, but leaves the small-scale structures intact. 

For some of our fits (see Sect. 12. 6K we use the spatially fil- 
tered quantities to define spatial masks, isolating certain mea- 
sured parameters in a given range. In particular, we use the fil- 
tered 538 nm continuum image to define masks that have intensi- 
ties in a specific range to estimate the relation between intensity 
and vertical and radial flows (see Sects. l3.ll and l3.2l) . We also use 
the filtered LOS magnetic field map to identify penumbral struc- 
ture where the LOS magnetic field is locally strong (spines) or 
weak and nearly horizontal (inter- spines). These masks, shown 
in Fig. [6l are based on the high-pass filtered ^los map, with 
thresholds set such that roughly 33% of the penumbra is con- 
tained in each of the two masks. 



2.6. Azimuthal fits 

The sunspot observed is at 15° heliocentric distance, such that 
the measured LOS velocity contains a mixture of contributions 
from both vertical and (strong) radial flowfl We must resort 
to statistical methods to separately determine average properties 
of the radial and vertical flows. The presently observed penum- 
bra has boundaries to the umbra and quiet Sun that are rather 
irregular, but the penumbral filaments and the radial extent of 
the penumbra appear similar on the disk center and limb sides. 
We assume that the small-scale velocity field and magnetic field 
fluctuations do not vary in amplitude with azimuth angle in order 
to separate the observed LOS component of the velocity into its 
radial and vertical components. 

Ignoring any azimuthal components of the flow field, the re- 
lation between the LOS velocity vlos, the radial (horizontal) 
velocity v^, the vertical velocity v^, the heliocentric distance 6 
(=15°) and the azimuthal angle is ther0 

vlos = COS (p sin ^ -h Vz cos 6 (4) 



(|Plaskettl[T952h . Fits based on this relation have been used to 
disentangle the radial and vertical components of penumbral 
velocity and/or magnetic fields fr om LOS measurements by nu- 



merous earlier investigators (e.g.,*Maltbv'l964":'Title et al. 


19931; 


SchHchenmaier & Schmidt l2000c .Bellot Rubio et al.l 


2003 


Tritschler et al.1 120041: iLanehans et al.1 12005; iBeUot Rubio et al.l 


20061 ISanchez Almeida et al.l 12007; IScharmer et al.ll201 1 


). The 



sign conventions here are that is zero in the disk center direc- 
tion, radial outflows are counted positive, and vertical velocities 
as well as velocities in the direction toward the observer are 
counted as negative. Our procedure for defining the azimuth 



^ A sunspot located precisely at disk center would be allow clean 
measurements of vertical velocities. However, even a few degrees helio- 
centric distance leads to significant contributions from radial Evershed 
flows that nearly alwa ys must be acco unted for in such measurements, 
as demonstrated bv aBharti et al.l (1201 Ih . 

^ The same relation, but with different sign conventions, is used to 
resolve 5los into its radial and vertical components. 



Fig. 6; Spine (upper panel) and inter- spine (lower panel) masks, 
obtained by thresholding the spatially filtered LOS magnetic 
field map, shown in the lower-right panel of Fig.[5j The thresh- 
olds have been set such that roughly 33% of the penumbra is 
contained in each of the two masks. 



angle and radial distances within the penumbra is described in 
SOM. The zero-point azimuth direction was determined from 
full disk images recorded by the Solar Dynamics Observatory 
(SDO) close to the time of our observations and is indicated by 
the arrows shown in, e.g.. Figs. [T] and [2l Guided by the orienta- 
tion of the penumbral filaments, we determined an approximate 
center of the sunspot from which straight lines line up as well 
as possible with the penumbral filaments. The umbra/penumbra 
boundary was identified by applying an intensity threshold 
and removing the light bridge and innermost part of a few 
filaments. From that boundary six radial zones were defined by 
successively applying a dilate operation using a circular kernel 
of 1.5'' radius. The boundaries of the six radial zones are shown 
as dark contours in Fig. \5\ Examples of azimuthal fits made 
at diff'erent high-pass filtered 538 nm continuum intensities 
can be found in Fig. [T] These show fundamental diff'erences 
between the spines and inter- spines: in the inter- spines the 
amplitude of the azimuthal variation of vlos implies a strong 
radial component of the LOS velocity that is much weaker in 
the spines and virtually absent at low continuum intensities. 
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Fig. 7: Examples of azimuthal fits of the absolute (unfiltered) 
LOS velocities to Eq. (4), measured with diff'erent methods in 
the 5380 and 6301 lines, for relatively dark (Sic ~ -0.19) and 
bright (Sic ~ 0.19) structures. The left column shows the data 
(dark dots) and fits (red curves) for spines, the right column for 
inter- spines. 
10 



Figures [3] and |4] show all measured quantities within two 5''x8'' 
subfields on (mostly) the limb and disk center sides, with tick 
marks separated at V intervals. The first and last panels in the 
top row show the 538 nm and 630 nm continuum images. These 
appear nearly identical but the continuum image at 538 nm 
shows somewhat higher spatial resolution, as expected. The first 
two panels in the bottom row show the 5380 COG and line 
core velocities. These are quite similar, but with significant dif- 
ferences for the smallest flow features shown. Although not so 
obvious in Fig. |4l there is an overall tendency, also outside the 
sunspot, for the COG velocity map to show stronger small-scale 
flows than for the 5380 line core velocity map. This is inter- 
preted to mean that the flows with the smallest scales change 
rapidly with height, whereas the large-scale flows show weaker 
gradients with height. 

The 6301 flow maps in general show much less fine structure 
than the 5380 flow maps. This is not a question of spatial reso- 
lution, as the continuum images demonstrate, but most likely is 
due to either or both of two efl'ects: the first is that small-scale up- 
flows do not extend to the same height as large-scale strong up- 
flows. The second effect is the relatively large formation height 
range of the 6301 line, aggravated by the limited spectral resolu- 
tion of CRISP, causing smearing of flow topologies that change 
with height within the formation range. We note that the 6301 
line core velocity (last panel) is the most "fuzzy" of these flow 
maps on the limb side, but not on the disk center side. We also 
note the 6301 COG and 70% bisector flow maps (third and fourth 
panels in the bottom row of Figs. [3] and |4]) appear almost iden- 
tical. In the following, we do not discuss these flow maps sepa- 
rately. We attribute this similarity to relatively low spectral res- 
olution of CRISP, smearing out the height variations of the flow 
field. 

We emphasize the small azimuthal scale of almost all the 
flow structures, in particular the (dark) red- shifted lanes in be- 
tween the (bright) blue- shifted features. Many of these flow 
structures appear to have azimuthal scales close to the SST 
diflfraction limit (0^.^14 at 538 nm). This clearly illustrates the 
difficulties of earlier (failed) attempts aimed at detecting con- 
vective downflows (i.e., that dark features are associated with 
downfl ows) in the interior penumbra at much lower spatial reso- 
lution (iBellot Rubio et aT]l2010l: iRanz & SchHchenmaied[2009l: 



lFranzll201 ll) . as demonstrated recently also from analysis of syn- 
thetic spectra calculated from simulation data and degraded to 
the spatial resolution of SST (Bharti et al. 2011). In the follow- 
ing, we characterize the complicated relations between the ob- 
served small-scale variations in the continuum intensity, flow 
field and magnetic field, and attempt to explain some of the ob- 
served differences between the limb and disk center sides of the 
penumbra. 

Finally, we note that the LOS magnetic field (mid panel, up- 
per row in Figs. [3] and |4]), shows only weak evidence for varia- 
tions that can easily be linked to the small-scale velocity field. 
However, the LOS magnetic field in Fig. |4l shows several nar- 
row dark structures, where either the polarity reverses or where 
the LOS magnetic field is nearly zero (according to the present 
COG measurements). These can in most cases be linked to lo- 
cally red-shifted ffow features. We will discuss these features in 
more detail in a forth-coming paper. 
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Fig. 8: Correlation between the high-pass filtered 538 nm continuum intensity and the high-pass filtered LOS velocity, measured with 
diff'erent methods in the 5380 and 6301 lines. The top row shows plots for the "interior penumbra", the bottom row corresponding 
plots for quiet Sun. For easier comparison, the red curves show in the upper panels are shown as dashed curves in the lower 
panels. Note that the velocities plotted are referenced to local averages by spatial filtering (unsharp masking). 



Table 1: Correlations and RMS velocities. 



Quantity 


5380 COG 


5380 core 


6301 COG 


6301 70% bis. 


6301 core 


R interior penumbra 


-0.686 


-0.610 


-0.693 


-0.602 


-0.404 


R quiet Sun 


-0.772 


-0.796 


-0.804 


-0.734 


-0.588 


Vrms interior penumbra (kms~^) 


1.08 


1.06 


0.58 


0.64 


0.32 


Vrms quiet Sun (kms~0 


1.38 


1.50 


0.91 


0.96 


0.56 



Notes. Summary of RMS velocities (vrms) and correlation coefficients {R) between high-pass spatially filtered 538 nm continuum intensities and 
LOS velocities (vlqs) measured with the 5380 and 6301 fines for quiet Sun and the interior penumbra. The different methods used for measuring 
the velocities are indicated above. 



3.1. Intensity and LOS velocity correlations 

A quick glance at Figs. [3] and |4] makes it evident that the mea- 
sured small-scale LOS velocities are strongly related to the con- 
tinuum intensity: e.g., blue-shifted (shown bright) features can 
in almost all cases be associated with a locally bright continuum 
feature. This does not necessarily mean that there exist simple 
one-to-one relations between these quantities: it is easy to find 
examples where a relatively strong blue- shifted feature shows 
only a marginal brightening in the continuum and vice- versa. In 
other cases, the flow features and continuum structures appear 
somewhat displaced spatially relative to each other. There is no 
reason to expect a perfect correlation between intensity (temper- 
ature) and flow velocities: even for quiet Sun granulation, this 
correlation is less than about 80%. We first quantify these cor- 
relations in the simplest possible manner. Figure [8] shows the 
correlation between the high-pass filtered 538 nm continuum in- 
tensity and the high-pass filtered LOS velocity for the interior 
penumbra (top row) and quiet Sun (bottom row). As shown in 
Fig. [21 the main eff'ect of high-pass filtering the LOS velocity 
map is to remove the systematic diff'erences in LOS velocity be- 
tween the disk center and limb sides - these diff'erences are from 
the radial (horizontal) flow. 

After high-pass filtering the LOS velocity map, the remain- 
ing small-scale flow field structures appear similar on the disk 



center and limb sides, suggesting a dominance of vertical flows 
(c.f., iBeckers & SchroteJl 19691 for similar arguments). The 
correlations between intensity and LOS velocity in Fig. [8] are 
obvious, the corresponding correlation coefficients are given in 
Table [Hand range from -60% to -69% (except in the 6301 line 
core, where it is -40%) for the penumbra. For the quiet Sun, the 
corresponding correlations are stronger and range from -73% 
to -80% (-59% at the 6301 line core). The values for the quiet 
Sun are similar to those found from bi sector measur ements close 
to the continuum in the 6301 line bv lFranzl (l201lh . -78%, the 
lower values found by him for the penumbra are likely to be 
explainable by the much lower spatial resolution of the Hinode 
data and the different methods of analysis used by Franz and us. 
Irrespective of these differences, our conclusion about the exis- 
tence of strong relations between the continuum intensity and the 
flow field does not rely on the precise values of the correlation 
coefficients found; as stated earlier they are obvious in Figs. [3] 
and[4l The correlation coefficients calculated here only quantify 
this relation in a simple way that is easily testable against sim- 
ulations. In the following we will analyze this in more detail, 
using an independent technique. We finally note (Table [B that 
the RMS velocities obtained for the penumbra are on the order 
of 65-70% of those for the quiet Sun, such that they suffice to 
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explain the penumbral heat flux (c.f. J Scharmer et all201 ll for a 
more detailed discussion). 

3.2. LOS velocity and magnetic field correlations 

Figure [9] shows the correlation between the measured LOS mag- 
netic field and LOS velocities in the interior penumbra, mea- 
sured with the COG method and from the line core of the 5380 
line, and the 70% bisector and from line core of the 6301 line. 
Note that the quantities shown in this figure have not been spa- 
tially filtered. All correlations shown (and that for the 6301 COG 
velocity) have a roughly triangular shape with the "corners" 
characterized as follows: The strongest blueshifts corresponds 
to weak (-400 G) LOS magnetic field with the same polarity as 
the leading polarity of the spot. The strongest redshifts also cor- 
respond to weak LOS magnetic fields but with a polarity that is 
opposite to that of the spot. Almost all this opposite polarity field 
is in the limb- side penumbra, but a few such patches can also be 
found in the center-side penumbra. Finally, the strongest LOS 
magnetic field (up to about 1600 G) systematically corresponds 
to small blue- and redshifts. 

In the lower-right panel of Fig. [TOl we show the variation 
of the radial velocity with the strength of the high-pass filtered 
LOS component of the magnetic field. This plot w as ob tained by 
creating LOS magnetic masks as explained in Sect. 12. 5 l an d m ak- 
ing azimuthal fits of the measured LOS velocities (Sect. 12.61) to 
separately determine the vertical and radial velocities. This plot 
shows a dramatic variation of the radial velocity with the filtered 
LOS magnetic field: where ^los is locally very weak, the radial 
velocity reaches up to 5 km s"^ at the height of formation of the 
5380 line, where it is very strong the radial velocity drops to less 
than 0.5 km s"M At heights corresponding to the formation of the 
6301 line, the radial velocities are systematically lower, but the 
overall trends are the same as for the 5380 line. The lower-left 
plot in Fig. [To] shows that the strength of the spatially filtered 
LOS magnetic field is closely related to the inclination of the 
magnetic field: where this quantity is small, the magnetic field 
is nearly horizontal and where it is strong, the magnetic field 
is inclined by about 50° to the vertical. Assuming that most of 
the variation in the (unfiltered) LOS magnetic field comes from 
inclination changes, rather than from field strength variations, 
these plots explain the triangular shape of the correlations shown 
in Fig.|9l Where the LOS magnetic field is weak (nearly horizon- 
tal field), strong radial flows produce considerable spread in the 
measured LOS velocities, with the extremes corresponding to 
the lower-left and lower right "corners" of the triangles. Where 
the LOS magnetic field is very strong, corresponding to more 
vertical field (upper "corner"), radial flows are virtually absent, 
such that the spread in the measured LOS velocities is limited to 
that of the (convective) vertical flows. 

In the remaining panels of Fig.[TOlwe show the results of fits 
of ^Los to the high-pass filtered ^los map. These plots show 
that selecting locally very weak LOS magnetic field corresponds 
to selecting a component of the penumbral field that is relatively 
weak (field strength of about 850 G) and (not surprisingly) al- 
most horizontal (inclination approximately 80°). Selecting the 
locally strongest LOS magnetic field corresponds to selecting 
the part of the magnetic field that is stronger (about 1200 G) and 
more vertical (inclination approximately 50° ). These propertie s 
agree with those of spines and inter-spines (iLites et al.lll993l) . 
although the much lower spatial resolution of their data com- 
pared to ours means that it is not obvious that we are iden- 
tifying the same penumbral structures. Nevertheless, we can 
establish that for the observed sunspot, the strong horizontal 
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Fig. 9: Correlation between LOS velocity, measured with diff'er- 
ent methods in the 5380 and 6301 lines, and the measured LOS 
magnetic field. Neither of the two quantities plotted have been 
high-pass filtered. 



(Evershed) flow prevails in the inter- spines and is virtually ab- 
sent in the spines. This division of the penumbra into two main 
components with quite difl'erent magnetic and flow properties 
is consistent with numerous results obtained from inversions 
based on spectrop olarimetric data at much lo \ yer spatial res- 
olution than here (ISolanki & MontayonI 1 19931: iMartmez Pillet 



20001: ISchlichenmaier & ColladosI |20Q1 iBel 



ot Rubio et al 



120031 120041: iBorrero et al.l 120071: iTritschler et al.ll2007h . and is 
at the center of interpretations in terms of a "flux tube" (at a 
smaller scale than our spines) and a "background component". 
The existence of these two fairly distinct populations is clearly 
illustrated in Fig.[TTl showing the 5380 COG and line core and 
the 6301 COG LOS velocity field within the inter- spine mask 
(top row) and within the spine-mask (bottom row). In the top 
row (inter- spines), we see clearly the signatures of the Evershed 
flow: strong blue- shifts on the disk center side and strong red- 
shifts on the limb side. In the spines, these signatures are missing 
entirely, showing that the Evershed flow is virtually absent. 

Figure [12] shows the radial variation of the average prop- 
erties determined for spines and inter- spines, indicating a 
magnetic field that dips down (inclination larger than 90°) 
in the inter- spines in the outer part of penumbra, consistent 
with earlier findings. In comparison, the plot of the average 
vertical velocity for the inter-spines shows a cross-over to 
positive velocities (downflows) further toward the center of 
the sunspot, clearly indicating that the average magnetic field 
and flow field are not precisely aligned (a similar result was 
obtained from the 5380 Stokes V and LOS velocities, see SOM 
Fig. SIO, lower-left panel). Systematic vertical and radial flows 
are small in the spines, but the inter- spines show strong radial 
flows up to 5 kms"^ in the interior penumbra and systematic 
upflows in the inner and downflows in the outer penumbra. 
Ignoring the existence of the small-scale inte nsity correlated 
vertic al (convective) flows demonstrated earlier (IScharmer et al.l 
I2OIII) and here, these flow patterns would be consistent 
with flows in arched flux tubes, as concluded in numerous 
earlier interpretations of observations made at much l ower 
spatia l resolution ( Solanki & MontayonI 1 19931: [M artinez Pilletl 
I2OOOI: ISchlichenmaier & ColladosI l2002l: iBellot Rubio etaLl 



12 



G.B. Scharmer and V.MJ. Henriques: SST/CRISP Observations of Convective Flows in a Sunspot Penumbra 




Fig. 11: Measured absolute (unfiltered) LOS velocities within the inter-spine (upper row) and spine masks (lower row) shown in 
Fig. [6l as obtained in the line core and with the COG method in the 5380 line and with the COG method in the 6301 line. The 
LOS velocities in the inter-spine mask show strong differences between the disk center and limb sides, demonstrating strong radial 
flows in the inter- spines, whereas no such difl'erences are seen in the spines. Tick marks are at V intervals. The arrow points in the 
direction of Sun center. 



2003L 120041: iBorrero et al.l 120071: iTritschler et aP l2007l: 
Franz & Schlichenmaieil l2009|) and numerical simulations 



(ISchlichenmaier et aL 1998; Schlichenmaier 20 02|). As pointed 
out earlier (ISanche z Almeida et al. 2007; Scha rmer et a 
and in Sect. 12. 2[ such apparent large-scale penumbral flow 
patterns are likely to be at least partly the result of inade- 
quate spatial resolution, giving efl'ects similar to convective 
blue- shifts. 



In Fig. [13] we show the results from azimuthal fits (Sect.[ 
of the measured LOS velocities made at diff'erent high-passed fil- 
tered 5380 continuum intensities, separately for spines and inter- 
spines. The fits include pixels only from the "interior penum- 
bra", defined here as radial zones 2-4 (see Figs. [T]-[2l and Fig.O. 
We find strong vertical velocities that are roughly proportional 
to Sic both in the spines and inter-spines, with the darkest struc- 
tures showing downflows of up to about 1 kms~^ (except in the 
line core of the 6301 line). The brightest structures show upflows 
up to about 3 kms"^ in the 5380 line and up to 1.5 kms"^ in 
the 6301 line. Strong radial flows are found only in the inter- 
spines. The 5380 COG measurements corroborate the earlier 
line c ore measurements made in the same line ( Scharm er et al.l 
l20TTh : here we detect similar clear convective signatures also in 
the 6301 line. The top four panels in Fig. [7] show considerable 



scatter around the fitted curves. This scatter shows only small 
variation with the azimuth angle and does not tend to disappear 
at = 90° and = 270°, where contributions from radial flows 
should be absent. This clearly indicates that the scatter comes 
from fluctuations primarily in the vertical flows that do not cor- 
relate with intensity. As concluded already from our data, the 
correlation between continuum intensity and vertical velocity is 
somewhat weaker in the penumbra t han in the qui et Sun. The 
same tendency is seen in simulations (lRempelll201lh . 

We make a few additional remarks relating to Fig. [131 First, 
the radial velocities obtained in the spines show a clear tendency 
to increase with brightness for the brightest structures. We note 
that Figs. [3] and in show several examples of bright strong flows 
that are located just at the boundaries of the spines. Our thresh- 
old for ^^L os set to d efine spines is "generous" compared to 
that used by lScharmer e t al. (201 1) based on the integrated 5380 
Stokes V data, where the spines occupied only 13.5% of the in- 
terior penumbra. Here, we set the threshold such that spines oc- 
cupy 33% of the penumbra. Increasing this threshold decreases 
the radial flow speed for the the brightest structures in the spines, 
and also reduces the strength of the radial flows at lower inten- 
sities. Second, the roughly linear relation between die and the 
average vertical velocity found will not change significantly if 
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Fig. 10: Variation of "interior penumbra" (radial zones 2-4) 
magnetic field and radial flow speed with the strength of the 
high-pass filtered LOS magnetic field, obtained from azimuthal 
fits. The top two plots show the variation of the vertical and hor- 
izontal components of the magnetic field, the lower-left the vari- 
ation in magnetic field inclination. The lower-right plot shows 
the (remarkable!) variation of the radial flow velocity, measured 
with five diff'erent methods from the 5380 and 6301 lines. Plus 
symbols correspond to measurements in the 5380 line, squares 
to the 6301 line, lines drawn full to COG measurements, short 
dashes to line core, and long dashes to 70% bisector (6301 line 
only) LOS velocities. The vertical dashed lines correspond to 
the upper and lower thresholds respectively, for "inter- spines" 
and "spines". 



the zero-point calibrations of LOS velocities are changed by one 
or two hundred m s"^ . This will change the value of the strongest 
downflows and upflows by the same amount but much larger er- 
rors would be needed to "remove" the penumbral dark down- 
flows. Such large errors would also imply that the penumbra has 
an average upward velocity of well o yer 0.5 kms ~^ in strong 
disagreement with earlier findings (c.f jFran3l201ll) . 

We finally note that the finding that the convective down- 
flows are of similar strength in the 5380 and 6301 lines, while 
the upflows are much stronger in the 5380 line than in the 6301 
line, suggests significant diff'erences in their height variations. 
This is consistent with the inversions of iRai^ (l20TTh . based on 
6301 line data obtained with Hinode from a sunspot close to disk 
center. He found that a typical upflow extend only up to heights 
corresponding to log(T) ^ -0.5 whereas typical downflows ex- 
tend up to heights corresponding to log(T) ^ -1.2. 




3 4 
Radial zone 



Fig. 12: Radial variation of properties of spines and inter-spines, 
as defined by the LOS magnetic masks, shown in Fig. (6] The 
plots in the top row show field strengths and inclinations for 
spines (full) and inter- spines (dashed). Symbols in the lower four 
plots are as in Fig.fTOl 



and the intermediate population). Except for the darkest struc- 
tures, the radial velocities are similar with and without stray light 
compensation and (as expected) intermediate to those of spines 
and inter-spines, shown in Fig. [131 The strongest vertical bright 
upflows are reduced in strength by only a small amount (from 
1.5 kms"^ to 1.3 kms"^) without stray light compensation, but 
the darkest downflows are reduced in strength by a factor of 2 
to 0.4-0.5 km s"^ . Given our estimated uncertainty of 150 m s"^ , 
we can thus establish the existence of the dark downflows from 
the 6301 line even without stray light compensation, but with re- 
duced margin for error in the zero-point calibration. We con- 
clude that both the high spatial resolution of the SST and com- 
pensation for straylight is needed for reasonably accurate esti- 
mates of the strengths of the dark convective downflows in the 
interior penumbra. 



3.3. Effects of straylight on measured vertical/radial flows 

Due to the strong weakening of the Ci 5380 line in dark gran- 
ular and penumbral structures, straylight compensation is cru- 
cial for detecting dark convective downflows in the penumbra 
in this line (see SOM). The efl'ects of straylight are less dra- 
matic for measurements of LOS velocities in the Fe i 6301 line. 
Figure [14] compares the vertical and radial velocities obtained in 
this line with and without straylight compensation. For these fits, 
we included the entire interior penumbra (spines, inter-spines 



4. Conclusions 

We have combined the measurements of penum bral LOS veloc- 
ities i n the 5380 line core reported on recently (IScharmer et al.l 
1201 ll) with LOS velocity and magnetic field measurements in 
the 6301 line as well as COG LOS velocity measurements in the 
5380 line. These observations were obtained with SST/CRISP 
at the highest spatial resolution (0^.^14 at 538 nm and 0^.^16 at 
630 nm) attainable today for spectropolarimetric data. By com- 
paring observed quiet Sun RMS continuum intensities and ve- 
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Fig. 13: Variation of vertical (top) and radial (bottom) veloci- 
ties with high-pass filter C i 5380 continuum intensity for spines 
(left) and inter- spines (right) in the interior penumbra (radial 
zones 2:4). Plus symbols correspond to measurements in the 
5380 line, squares to the 6301 line, lines drawn full to COG 
measurements, short dashes to line core, and long dashes to 70% 
bisector (6301 line only) LOS velocities. 
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Fig. 14: Comparison of vertical (top) and radial (bottom) veloc- 
ities measured in the 6301 line, with (left column) and without 
(right column) stray light compensation. All structure in the in- 
terior penumbra (spines, inter-spines and the intermediate pop- 
ulation) were included in the fits. Lines drawn full correspond 
to COG measurements, short dashes to the line core, and long 
dashes to the 70% bisector LOS velocities. 



locities to values established from simulations of convection, we 
estimate the level of straylight in our data. Constraints on the 
FWHM of the straylight PSF are obtained from the minimum 
umbral intensity, which is low (15-18%) in the observed data. 



By using these constraints, we compensate the data for spatial 
stray fight (IScharmer et al .1120 111) . 

We assume that the measured small-scale fluctuations in in- 
tensity and LOS velocity across penumbral filaments do not 
vary systematically with azimuth angle at constant radial dis- 
tance from the inner boundary of the penumbra. Azimuthal 
fits of the measured LOS velocities can then be used to deter- 
mine statistical properties of vertical and radial flows separately 
( Plaskett I952I ). We have extended this widely used method (e.g. , 
Maltbv 1964; Title et al. 1993; Schlichenmaier & Schmidt 2000'; 

Tritschler et al 2004: Langhans et al. 2005: Bellot Rubio et al] 

[2iD06l: Sanchez Almeida et al. by making azimuthal fits 

separate ly for diff'erent (s patially filtered) continuum intensi- 
ties ( Sch armer et al. l l20Tll) . These fits are made for the inte- 
rior penumbra (radial zones 2-4), excluding its innermost and 
two outermost radial zones (note that th e previous analysis in - 
cluded radial zones 2-5, c.f. SOM in Scharmer et al. r i2oTii) . 
We note that 'Sanchez Almeida et al. I (l2007l) made azimuthal fits 
by assuming linear relations between local intensity fluctuations 
and vertical/horizontal velocities. These fits demonstrated a clear 
correlation between intensity and the vertical velocities in the 
penumbra, in the sense expected for convection. However, this 
method cannot be used to draw conclusions about whether or 
not the dark component of the penumbra on the average is asso- 
ciated with actual downflows (v^ > 0). 

Our recent analysis of line core measurements in the 5380 
line demonstrated a roughly linear relation between the penum- 
bral continuum intensity and the vertical velocity (SOM, Fig. 
SIO). We also found that the darkest structures on the average 
show downflows of about 1 kms~^, whereas the brightest struc- 
tures show upflows of up to 3 km s"^ in t he layers correspondin g 
to the formation height of the 5380 line dScharmer et al.ll201ll) . 
Here, we find similar downflow velocities in the dark struc- 
ture but only about 1.5 kms"^ bright upflows in the wings of 
the 6301 line. This analysis corroborates the results obtained, 
and conclusions draw n from the 5380 line core measurements 
(IScharmer et al.ll20TTI) . but indicates that the strong bright up- 
flows seen in the 5380 line do not extend to the heights of the 
6301 line. We also find that the convective signature for vertical 
flows is the same irrespective of whether the magnetic field is 
relatively vertical and strong (in the spines) or weaker and more 
horizontal (in the inter- spines). However, strong radial outflows 
are found only in the inter- spines. This suggests the existence of 
either two distinct modes of (magneto-) convection simultane- 
ously in the penumbra, where the strong Evershed flow is related 
to convection only in the more horizontal and weaker magnetic 
field structures in the inter- spines, or a more gradual change of 
the horizontal flow topology with magnetic field inclination. 

Correlations between local fluctuations in intensity and LOS 
velocity in the penumbra in the sense expected for convec- 



tion have been reported repeatedly in the literature since 1969 
(e.g.. iBeckers & Schroted 119691: ISanchez Almeida et al." '19931: 


Schfich 


lenmaier & Schmidt 


1 19991: ISchmidt & Schlichenmaier 


20001; 


Sanchez Almeida et al.l 120071: iFranz & Schlichenmaier 


|2009|; 


Franzl I2OIII). However, the interpretation of 


these 



observations in terms of convective flows is not without 
ambiguity. At low spatial resolution, the large-scale struc- 
ture of these flows appears to be that of an upflow of 
bright structure in in the inner penumbra and a down- 
flow of dark structure in the outer parts. This flow pat- 
tern has been interpreted as support for flux tube models 
( Schlichenmaier & Schmidt 1999; Sch midt & Schlichenmaieri 
.2000: .Franz & SchlichenmaieB2009:.Franzll201ll) . 
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A major objection to the interpretation of these flows 
as being of convective origin is the absence of evidence 
of dark downflows in the in t erior body of the penumbra 
(iFranz & Schhchenmaied l2009l: iFranzl l201lh . In our analysis, 
we therefore exclude the outermost penumbra, where dark 
downflows have been reported repeatedly, from our analy- 
sis. We also exclude the innermost part of the penumbra, 
where bright upflows are omnipresent and have been inter- 
preted as the footpoints of penumbral flux tubes. The discov- 
ery of d ark (convective) downflows in the penum bra reported 
recently (IScharmer et al.l l201 1[ fjoshi et al.l l201 ih and the cor- 
relation bet ween intensity and vertical flows found from the 
5380 line (iScharmer et al.l 120111) in the interior penumbra, is 
confirmed here for the 6301 line and for COG velocities in 
the 5380 line. The presence of intensity correlated upflows and 
downflows at roughly the same radial distance in the interior 
penumbra contradicts earlier interpretations in terms of embed- 
ded flux tubes. The present observations, combined with earlier 
past and recent indications from correlatio ns between penum- 
bral brightness and LOS velocities (e.g., iBeckers & SchroteJ 
119691: ISanchez Almeida et al.l l2007b. and observed dynamics 
of penumbral filaments (llchimoto et al.l l2007bt IZakharov et al.l 
120081: ISpruit et aDl2010l) . provides overwhelming observational 
evidence for interpreting the filamentary structure, dynamics, 
heat flux, magnetic field topology and Evershed flow of the 
penumbra as due to convection, rather than in terms of em- 
bedded flux tubes . This interpretation is consistent with theo- 
retica l arguments (ISpruit & Scharmei] |2006^, 'Scharmer & Spruitl 
2Q06h as well as numerical simulations (Heinemann et al. 2007; 
Scharmer et ani2008bl: iRempel et al1l20()9aiibl : IRempelll2011 ). 

Our observati ons could also be i nterpreted as evidence for 
convection rolls dPanielsonl [l961bl). as has been suggested 
previously (i Beckers & Schroter* '1969'; 'ich imoto et al ] I2007bl: 
IZakharov eF al. 2008; Sanchez Almeida 201 3) -However, objec- 
tions against this interpretation were raised by lScharme3 (l2009l) 
on the following grounds: the convection rolls are expected to 
be associated with a n early horizontal magne tic field, whereas 
e.g., the observations of llchimoto et al.l (l2007bl) refer to the inner 
penumbra, where the magnetic field has a strong vertical com- 
ponent. Also, to sustain the radiative output of a filanient oy er its 
life time of on the order of an hour (iLanghans et al.ll2QQ7b , the 
convective upflows must persist to depths much larger than a few 
100 km, whereas the convective roll is a shallow phenomenon 
d^nielson 1961b). 

The combination of high spatial resolution and straylight 
compensation of our data results in RMS velocities for the 
penumbra that are on the order of 65-70% of those for the quiet 
Sun, such that th ey suffice to explain the penumbral heat flux 
(c.f. JScharmer et al. 201 1, for a more detailed discussion). This 
also strongly suggests that the fundamental scales at which this 
convection occurs is resolved in our data. 

Our data also contain a wealth of information about the 
toplogy of the convective flows and their magnetic fields. In 
future work, we will investigate this in more detail. We will 
also compare the present data with simultaneously obtained even 
more highly resolved images in the Ca H line (Henriques, in 
prep.), with emphasis on exploring properties of bright penum- 
bral grains. We also intend to investigate properties of penum- 
bral dark cores (iScharmer et al.t2002 ), which are obvious mostly 
in the innermost penumbra (radial zone 1), and thus were to a 
large extent excluded in the present analysis. We do not exclude 
the possibility, however, that a small fraction of pixels within 
our mask defining the interior penumbra actually corresponds to 
dark cores, and thus are included in our azimuthal fits. 
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